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A general scenario that leads to Coulomb quantum crit- 
icality with the dynamical critical exponent z = 1 is pro- 
posed. I point out that the long-range Coulomb interaction 
and quenched disorder have competing effects on z, and that 
balance between the two may lead to charged quantum crit- 
ical points at which 2 = 1 exactly. This is illustrated with 
the calculation for the Josephson junction array Hamiltonian 
in dimensions D = 3 — e. Precisely in D = 3, however, the 
above simple result breaks down, and z > 1. Relation to other 
studies is discussed. 

The crucial difference between the quantum (T = 0) 
and the more familiar finite temperature phase transi- 
tions is that while dynamics is irrelevant for the latter, 
it is essential for the former The link between stat- 
ics and dynamics at a continuous quantum phase transi- 
tion is usually parameterized with the value of dynamical 
critical exponent z that describes relative scaling of the 
time and the length scales in the problem 1^, 1^. To- 
gether with the correlation length exponent v, z enters 
the low-temperature scaling of all physical observables, 
since in the vicinity of a quantum critical point temper- 
ature scales as r ^ I'^l^", where 8 is the T = tuning 
parameter. The value of z for a given quantum critical 
point is therefore of great interest, and has often been 
used to distinguish one universality class from another. 
In their seminal paper on universal conductivity in two 
dimensions. Fisher, Grinstein and Girvin Q also pro- 
posed that when the long-range Coulomb interaction is 
present, at the criticality energy should scale as inverse 
of length, so that z — \ should result. This well-known 
conjecture has since been used in interpreting some of the 
most intriguing experiments in modern condensed mat- 
ter physics, ranging from superconductor-insulator tran- 
sitions in low-dimensional systems Q , via metal-insulator 
transitions in Si-MOSFET's 0, to the universality of the 
underdoped high-Tc cuprates^ . It has also been utilized 
in the Monte Carlo simulations where knowing z in ad- 
vance greatly simplifies the inevitable finite-size scaling 
analysis 

The purpose of this Letter is to provide the theoretical 
justification for this widely used relation and point to its 
limitations. I show that in certain dimensions Coulomb 
coupling constant (i. e. charge) is protected from renor- 
malization, and consequently its flow under scaling trans- 
formation is directly proportional to its canonical dimen- 
sion, which is just z—1 [^1 (see Eq. (2)). This implies that 
i/ a charged critical points exists in the theory, its loca- 
tion is actually determined by the solution of the equation 



z = 1, which then also determines the value of dynami- 
cal critical exponent exactly. I argue that this situation 
arises when there is an additional coupling in the theory 
with the competing effect on the dynamical exponent, 
and which can balance the effect of Coulomb interaction. 
Such a coupling is shown to be provided by quenched 
disorder, and a concrete realization of the above scenario 
is worked out on the example of a disordered Josephson 
junction array Hamiltonian in Z? = 3 — e dimensions. Fi- 
nally, the simple relation z = 1 is found to break down 
at special dimensions at which the above renormalization 
group (RG) protectorate on charge is lifted. Relation to 
other recent theoretical studies of Coulomb criticality is 
discussed. 

To be specific, I will focus on the theory originally 
considered in [|l^ which describes an array of coupled 
Josephson junctions, but it will transpire that the un- 
derlying physics is more general. Building on an earlier 
work by Ma Fisher and Grinstein |jl^ have shown 
that Coulomb interaction can be represented by a min- 
imal coupling to the soft scalar gauge field. In the long 
wavelength limit the critical field theory that describes 
the system of bosons interacting via Coulomb interaction 
at T = and in /^-dimensions takes the form 
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where by |V|^^^ it is meant |fc|^^^ in the Fourier space, 
^'(r, r) is the complex superfluid order parameter, and 
Ao{f,T) is the scalar gauge field, which when integrated 
out introduces the Coulomb interaction into the remain- 
ing action for 4". I have also included a random poten- 
tial V{f,T), more on which in a moment. More gener- 
ally, the gauge-field propagator in the momentum space 
is (e^(K(fc) — 1)) @, which for the Coulomb interac- 
tion Vc(r) = l/r then yields the last term in (1) at small 
momenta. The presence of a neutralizing background is 
included by omitting the fc = components of the gauge 
field §. 

The above action without randomness (V(x, t) = 0) 
was first studied by Fisher and Grinstein [|lO| , and re- 
cently revisited by Ye jl^. In D = 3 both the charge 
(e^) and the quartic interaction (A) are marginally irrele- 
vant, and the Coulomb interaction causes the fiow to run 
away to negative A, which may be interpreted as a sign 
of a discontinuous transition. In D = 3 — e the result 
is more interesting: charge is irrelevant at the XY crit- 
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ical point, but if too large may still lead to a runaway 
flow. Irrelevance of the charge in the theory (1) arises 
in a somewhat non-trivial way, and it will prove instruc- 
tive to understand this in some detail. In Z? < 3 the 
inverse gauge-field propagator in (1) is non-analytic in k 
and therefore can not get renormalized by the integra- 
tion over the high-energy modes. This RG protectorate is 
reminiscent of the situation in the 2 + 1-dimensional elec- 
trodynamics with the Chern Simons term , where the 
statistical angle is exactly marginal for a similar reason. 
Here, this means that the /3-function for the charge is de- 
termined exclusively by its canonical dimension. Gauge 
invariance and the accompanying Ward identity imply 
that Aq ^ t^-^, so assuming t ^ yields to 
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in Z) < 3, where b is the standard RG parameter. The ex- 
ponent z in the last equation then needs to be determined 
from the renormalization of the ^'-propagator. The sim- 
plest one-loop calculation (see the first diagram on Fig. 
la))inr' = 3 — e then gives z = 1 — e^/3, and thus 
smaller than one. By the Eq. (2) small charge is then 
irrelevant. Two features of this result that are likely to 
be quite general should be noticed: first, z 1 since the 
scalar gauge-field couples only to the time derivative, and 
thus discriminates between space and time. Secondly, the 
negative sign in the Eq. (2) comes from the tendency of 
the gauge-field to soften the ^'-propagator. It is analo- 
gous to the well-known (and much debated) result that 
the anomalous dimension in the scalar (Higgs) electro- 
dynamics is negative ||l^, [Q, jl^. It could therefore 
be expected that when higher order terms in are in- 
cluded z < 1 will remain, and that irrelevancy of the 
charge in the pure theory (1) is more general than the 
simple one- loop calculation would suggest. 

The Eq. (2) also implies that if there would exist 
a critical point in the theory with non-zero charge, in 
D < 3 a,t that fixed point z = 1 exactly. This non- 
perturbative result parallels another exact result i]a = ^ 
in the D = 3 scalar electrodynamics , [|l6| where tja is 
the anomalous dimension of the vector gauge-field. The 
crucial question is how can such a Coulomb critical point 
(with ^ 0) arise. A clue is provided already in the 
above discussion: one needs another coupling which will 
tend to increase z and balance the effect of Coulomb in- 
teraction in the Eq. (2). A physically realistic candidate 
is disorder: in the quantum theory the random potential 
y(r, r) is random in space but static and independent 
of imaginary time. This anisotropy will in general lead 
to a non-trivial z, and it is well-known that the effect to 
the lowest order is always to increase it |p^ . In other 
words, while the weak Coulomb interaction is marginally 
irrelevant at the XY critical point, I expect it to become 
relevant at the disordered critical point where z > I. 
Small charge near the random critical point should grow 
until it balances disorder in the Eq. (2). As a result, a 



(a) 



(b) 















+ 





























(c) 



"0~ ^ Q 



(d) 



FIG. 1. One loop contributions to the renormalization of 
self-energy, quartic interaction, disorder coupling, and polar- 
ization. Wavy lines represent the gauge-field propagator, and 
dashed lines the disorder vertex. 

new stable (Coulomb) critical point may arise, at which 
z = 1 exactly. 

Next I demonstrate that the above scenario is indeed 
born out in the one- loop calculation in the theory (1), in 
passing reconciling the apparently different results in refs. 
pOf and Q. I then proceed to show how the equality 
z = 1 is violated in D = 3, and comment on relations to 
other works. 

To exert some control over the fixed points in the the- 
ory I will assume both a small e ~ 3 — D, and e^, where 
the latter is the number of dimensions over which dis- 
order is correlated The physically interesting case 
correspond to e^. = 1, but since I am primarily inter- 
ested in the point of principle, convergence properties 
of the double-e expansion will not be of much concern 
here . To average over disorder I utilize the standard 
replica trick, which attaches a replica index a = 1, ...TV 
onto all fields in the action (1) and introduces another 
interaction-like term in the theory: 
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is to be taken at the end of 
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the calculation ||l8[. Here R = (f, r), and for = 1 
one recovers the quantum problem where disorder is un- 
correlated (Gaussian) in space and independent of the 
imaginary time. 

To perform the Wilson-Fisher momentum shell renor- 
malization group one integrates out all the fields with 
A/6 < A: < A and — cxo < uj < oo, where A is the ul- 
traviolet cutoff. The effect of this procedure is to alter 
the coefficients in front of tj^ (Z^), fc^ (Zk), rn? A 
{Zx), W {Zw), and ^Q-term (Za) in the Fourier trans- 
formed action (1) |^^. One then rescales the momenta 
and the frequencies as 5fc — *■ /c and b^u — > to, and the 
fields as b-'^Ao Aq and b~ ° ^ Zl^^^/ to 
find finally that by defining new coupling constants as 
A(6) = b'^-^-'Z-^ZxX, W{b) - b^~^-'+'^Z-^ZwW, 
and e^(6) = b^~^Z^^e^ the action can be restored into 
its original form (at the critical point = 0). Com- 
puting next the Z-factors diagrammatically to one-loop 
order (Fig. 1) gives at the criticality 
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The exponent z in the Eq. (4) is determined by demand- 
ing that b^'^^'Z^ = b^^Zk, which gives 



z = l 
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Note that for D < 3 the fiow equation for the charge 
reduces to the Eq. (2). Precisely in D = 3 the inverse 
propagator for the gauge field becomes analytic, ^ . 
This means that in Z? = 3 the charge becomes renormal- 
ized by the polarization diagrams in Fig. 1 d), which 
contribute the last term in the Eq. (4). For W = the 
above /3- functions reduce to those of ref. (l^ when D < 3, 
and coincide with those of ||l^] right at D = 3, upon sim- 
ple redefinitions of couplings. They are also equivalent to 
those of for = 0. It may be interesting to note that 
many of the individual diagrams on Fig. 1 are ultravio- 
let divergent, due to the independence of the gauge-field 
propagator on frequency. All those divergences exactly 
cancel out in the final result Q. Finally, the flow of the 
mass term in (1) yields the correlation length exponent: 



1 1,, e2 W. 
u = - + -(A H ). 

2 4^6 4 ' 
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Let us turn now to the fixed points of the above equa- 
tions. Besides the Gaussian and the XY fixed points 
at W = = 0, both unstable with respect to dis- 
order, there are two disordered fixed points. First, at 



= 0, there is a neutral disordered fixed point at 
A„ = 2(4e + ner)/9 and Wn = 8(e + 5e^)/9, which is 
attractive in the X — W plane. At the neutral fixed point 
Zn = 1 + W/8 > 1, SO a weak Coulomb interaction is a 
relevant perturbation. With a small charge turned on, 
the fiow is towards a new, stable. Coulomb critical point. 
In Z? < 3 this fixed point is located at 
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Wc = —i2X,~e~er), 



and 
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where the last equation ensures that z = 1. The reader 
should note that disorder is necessary for the existence of 
the Coulomb critical point: without it the critical point 
would turn imaginary, and one would find only the stan- 
dard runaway fiow characteristic of the gauge-field fiuc- 
tuations in the e-expansion [|l[. The Coulomb critical 
point therefore may be considered as an example of a 
disorder induced continuous phase transition. 

Precisely in Z? = 3 the last term in the Eq. (4) becomes 
finite. There still exists a stable Coulomb critical point 
at Xc = 3.6Ur, Wc = 40(2Ac-e^)/41, and el = 3M^c/10, 
but with the dynamical critical exponent 



(12) 



which gives z « 1.15, for example, for Ct = 1. One finds 
2;7^1in£' = 3asa result of the removal of the RG pro- 
tectorate on charge. The same violation of the simple re- 
lation z = 1 can be expected in other problems in special 
dimensions. Also, since at the criticality Vc{r) ^ ^/r^, 
one expects that z > 1 in general [Q, since screening 
should certainly not make the interaction longer ranged. 
The result (12) thus implies the Coulomb interaction has 
been partially screened at the criticality in Z? = 3, and 
now decays faster (but still as a power-law) with distance 

H- 

The result z = 1 was previously also found in the 
large-N theory of Dirac fermions in D — 2 j2^, inter- 
acting both via Coulomb interaction and with the Chern 
Simons field. This theory may be relevant to the quan- 
tum Hall state-insulator transition ||^ , |2j] . It was found 
that for a certain range of the statistical angle (Chern 
Simons coupling) the competition between the Chern Si- 
mons and the Coulomb interaction leads to a non-trivial 
charged fixed point, at which z = 1 exactly, just like in 
the Eq. (2). The physical reason why the two interac- 
tions have competing effects on z only for some values 
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of the Chern Simons coupling remained somewhat ob- 
scure in this work. Nevertheless, the result in |^ bears 
a formal resemblance to mine. 

Finally, z = 1 was also found in the previous work 
by the author on the quantum critical behavior of dirty 
bosons with Coulomb interaction in 13 = 1 + e dimensions 
psf . There it arises as a consequence of a special symme- 
try the theory dual to (1) possesses precisely in D = 1, 
and can therefore be suspected to be an artifact of the 
specific RG scheme that was employed. Present paper 
can thus be understood as complementing the previous 
work in that it shows that z = 1 is also exact near the 
D = 3, and it may therefore be expected to hold in the 
physical case D = 2. 

It may also be interesting to note that simply setting 
e = er = 1 in the one- loop Eqs. (7) and (8), to crudely 
estimate the exponents in Z) = 2, besides the exact z ^ 1 
also yields v = 1.46 at the Coulomb criticality. Exper- 
imentally, « 1 jsj, and the result for v is less reliable 
than the one obtained in the expansion around D = \ 

H- 

To summarize, it was shown that the competition be- 
tween the correlated (quantum) disorder and Coulomb 
interaction may lead to new charged critical point at 
which the dynamical critical exponent z = 1 exactly. 
This simple result breaks down in special dimensions, 
and an example of the Josephson junction array when 
this happens in 13 = 3 was provided. I argued that sim- 
ilar results should be expected whenever there are two 
couplings in the theory with competing effects on z. 

The author is grateful to Matthew Fisher and Shiv- 
aji Sondhi for their useful comments and interest in this 
work. This research was supported by NSERC of Canada 
and the Research Corporation. 
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